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Over the last ten years tri-bimaximal mixing has played an important role in modeling the flavour problem. 
We give a short review of the status of flavour symmetry models of neutrino mixing. We concentrate on 
non-Abelian discrete symmetries, which provide a simple way to account for the TBM pattern. We discuss 
phenomenological implications such as neutrinoless double beta decay, lepton flavour violation as well as 
theoretical aspects such as the possibility to explain quarks and leptons within a common framework, such 
as grand unified models. 
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1 Introduction 

The long-standing solar and atmospheric neutrino anomalies suggested the idea of neutrino oscillations, 
now confirmed by a series of "laboratory" experiments based on reactors and accelerators [ 1 1. Especially 
puzzling theoretically is the fact that the neutrino mixing angles inferred from experiment follow a pattern 
rather different from that which characterizes quark mixing |2,3|. The atmospheric angle 823 is close 
to maximal, with a large value for the solar angle #12, both of which are at odds with their quark sector 
counterparts. Moreover, following the first indications of nonzero 6*13 reported by accelerator experiments 
MINOS [4J and T2K |5j three recent measurements of 6* 13 have been reported by the reactor experiments 
Double CHOOZ Daya Bay Q and RENO |§], as well as by the MINOS collaboration g [] 

While the historic discovery of neutrino oscillations provides strong indications for the need of physics 
beyond Standard Model (SM), the detailed nature of this physics remains elusive: i) the mechanism re- 
sponsible for neutrino mass generation, ii) its flavour structure, iii) its characteristic scale, as well as the 
nature of the associated messenger particle all remain unknown. As a result the nature of neutrinos, their 
mass and mixing parameters are so far unpredicted flO) . 

Understanding the pattern of neutrino mixing is part of the flavour problem, one of the deepest in particle 
physics. Although it may be the result of an accident, the pattern of neutrino mixing angles most likely 
follows a rationale. Indeed there has been a strong effort towards the formulation of symmetry-based 
approaches to address the flavour problem from first principles, assuming the existence of an underlying 
flavour symmetry of leptons and/or quarks, separately or jointly. 

In 2002 Harrison, Perkins and Scott proposed the tri-bimaximal (TBM) mixing ansatz fTT) with ef- 
fective bimaximal mixing of and v T at the atmospheric scale and effective trimaximal mixing for v e 
with and v T at the solar scale (hence 'tri-bimaximal' mixing). While large atmospheric mixing was 
already discussed before 2002, the trimaximal solar angle has represented a milestone for model building. 
Non-Abelian continuous and discrete flavour symmetries have been extensively used to account for TBM 
mixing. Here we review the basic features of some of the most interesting models proposed in the last ten 
years. 

To be fair we must say that the global analysis of neutrino oscillation data now indicates a robust 
measurement of a relatively "large" value of 6*13 JT5} which casts some doubt on the validity of the TBM 
ansatz as a good first approximation to the neutrino mixing pattern. Nevertheless it is too early to jump 
into conclusions, since in concrete theories there may be large corrections to the TBM pattern, so that here 
we still take it as a useful reference ansatz- 

Non-Abelian discrete groups have non trivial irreducible representations (irreps). Assigning the three 
known generations of leptons to irreps of a flavour symmetry group one can make predictions for masses 
and mixings in the lepton sector. In general it is expected that the number of free parameters of models 
based on Abelian flavour symmetries is typically larger then the corresponding number of free parameters 
needed to describe non- Abelian flavour symmetry models. Moreover there are non- Abelian discrete groups 
that contain triplet irreps, exactly as the number of generations in the standard model. Hence there are 
viable and predictive non- Abelian models to which we dedicate this brief review. 

The smallest group^Jthat contains triplet irreps is A4, the group of the even permutations of four objects, 
isomorphic to the group of the symmetries of the tetrahedron T. For a classification of the irreps of 
different non Abelian discrete groups see for instance fPJ) . A4 was first used in the lepton sector by 
Ma and Rajasekaran [14| but the solar angle was not predicted and neutrino masses were degenerate. A 
realistic model was proposed by Babu, Ma and Valle [15| adopting a supersymmetric context in order to 
produce the required neutrino mass splittings and mixing angles, predicting maximal atmospheric mixing 
and vanishing #13 to first approximation. Although expected to be sizeable, the solar mixing angle is 
unpredicted. In order to predict the full tri-bimaximal pattern, the neutrino mass matrix must take the form 



The bulk of the data on neutrino oscillations are well described in terms of three active neutrinos. 
The order of a finite group is just the number of elements. 
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I y x x \ 

M v = \xy + zx — z\, (1) 

\ x y — z x + z J 

where x, y, z are free parameters. The above matrix has two properties: 

• it is fi — t invariant giving maximal atmospheric and zero reactor angles; 

• it satisfies the relation (M„)n + (M„)i2 = (Mv)ti + {Mv)n giving trimaximal solar agle. 
The neutrino mass matrix of eq. ([TJ is diagonalized by the TBM mixing matrix 

/ 2/^6 1/V3 \ 
U= -1/V6 l/v/3 1/72 , (2) 
V -1/V6 1/V3 -l/x/2 / 

independently of the mass eigenvalues. A trimaximal solar angle was first given in a paper of Ma fl6) 
based on type-II seesaw but assuming in an ad hoc way that the contribution of two scalar singlets 1 and 
1' were the same. The derivation of TBM mixing from a flavour symmetry was achieved by Altarelli and 
Feruglio in Refs. fI7) and |T8) and subsequently by Babu and He fl9) . 

In contrast to the quark sector, neutrino mixing angles are large, possibly TBM, since the flavour group 
can break into two different subgroups in the charged and neutral lepton sectors, respectively. Consider 
A4 as an example. A 4 contains two abelian subgroups, namely Z 2 and Z 3 . When broken into Z 3 in the 
charged sector, and into Z2 in the neutrino sector, A4 leads to a lepton mixing matrix of TBM form, Eq. (|2}. 

Of course A4 is totally broken, therefore deviations at next to leading order are expected. In general 
one can not align A4 in the Z3 and Z2 directions in the charged and neutral lepton sectors respectively, 
this is known as the alignment problem. This may be circumvented by using extra dimensions and/or 
supersymmetry fT7|[T&) or by assuming a suitably chosen soft breaking sector. Alternatively, using a large 
discrete group, namely Z\ x {/^(l) 3 x S3, Grimus and Lavoura have shown |20| how to obtain the TBM 
form without alignment problem. 

2 The origin of neutrino mass 

Table [T] lists the fifteen fundamental "left-handed" chiral fermions of the Standard Model (SM), sequen- 
tially repeated, one set for each generation. In contrast to charged fermions, neutrinos come only in one 
chiral species, and parity violation in the weak interaction is introduced explicitly by having only "left" 
fermions transforming as doublets under the SU (3) ® SU(2) ® U (1) gauge group. 

The simplest and most general way to generate neutrino mass in the Standard SU(3) <g> SU(2) ® U (1) 
Model (SM) is by adding an effective dimension-five operator O a b = ^abf-a^b where l a denotes any of 
the three lepton doublets and $ is the SM scalar doublet. pT) . After electroweak symmetry breaking takes 
place, through the nonzero vacuum expectation value (vev) (<£>), Majorana neutrino masses are induced. 
From such general point of view the emergence of Dirac neutrinos would be an "accident", justified only 
in the presence of a fundamental lepton number symmetry, in general absent. The underlying nature of the 
dimension five operator in Fig.[T]is unknown: little can be said from first principles about the mechanism 
that engenders O a b, its associated mass scaJe or its flavour structure. The strength of the operator O a b can 
be naturally suppressed if the associated messengers are superheavy, as expected say, in unified scenarios. 
Alternatively, its strength can be naturally suppressed even in the absence of heavy messengers, due to 
the fact that O a b violates lepton number by two units (AL = 2), i.e. in its absence the theory recovers 
lepton number conservation. This is known as t'Hofft's naturalness [22 1 . Correspondingly, one may have 
high and low-scale neutrino mass models, depending on the mass characterizing the messengers whose 
exchange induces O a b- While the former type are closer to the idea of unification, the latter are closer to 
experimental testability. 
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Fig. 1 Dimension five operator yielding neutrino mass. 
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Table 1 Lepton, quark and scalar multiplets of the Standard Model 



2.1 High scale seesaw mechanisms 

The exchange of heavy messenger states, either fermions (type-I or type-Ill seesaw) or scalars (type-II 
seesaw) provides a simple way to generate the operator O a b- The smallness of its strength is ascribed to 
the large mass scale characterizing the violation of total lepton number |23|p4| . The simplest and most 
general description of the seesaw mechanism is in terms of just the SU (3) <g> SU(2) <g> U(l) gauge group 
with ungauged lepton number broken either explicitly J25| or spontaneously p6| . The latter framework 
or "1-2-3" scheme is characterized by SU(3) (g> SU(2) <g> U{1) singlet, doublet and triplet mass terms, 
described by the matrix [25 , 26) 

M,= (^l I"" 2 ) 0) 

where v 2 = ($) denotes the SM Higgs doublet vev and the basis is v^, z/£, corresponding to the three 
"left" and three "right" neutrinos, respectively. Note that, though symmetric, by the Pauli principle, M. v is 
complex, so that its Yukawa coupling sub-matrices Y v as well as Y3 and Yj are complex matrices denoting 
the relevant Yukawa couplings, the last two symmetric. Such SU(S) ® SU(2) ® U(l) seesaw contains 
singlet, doublet and triplet scalar multiplets, obeying a simple "1-2-3" vev-seesaw relation of the type 

v 3 vi ~ v 2 2 with Vx > v 2 > v 3 (4) 

This vev-seesaw is consistent with the minimization condition of the SU(3) <8> SU(2) ® U(l) invariant 
scalar potential, and implies that the triplet vev V3 — > as the singlet vev v\ grows. Neutrino masses are 
suppressed either by heavy 577(3) ® SU(2) ® U (1) singlet "right-handed" neutrino exchange (type I) or 
by the smallness of the induced triplet vev that follows from heavy scalar exchange (type II), as illustrated 
in Fig. [2] The matrix M.^ is diagonalized by a unitary mixing matrix U v , 

UlMvUv = diag(mi, M;), (5) 

yielding 6 mass eigenstates: the three light neutrinos with masses rrii, and the three heavy two-component 
leptons. The light neutrino mass states V{ are given in terms of the flavour eigenstates via the unitary matrix 
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Fig. 2 Left: type-I seesaw (right-handed neutrino exchange). Right: type-II seesaw (scalar triplet exchange). 



(6) 



o=l 



where the diagonalization matrices are given as a perturbation series, see Ref. |26|. The effective light 
neutrino mass, obtained this way is of the form 



Y 3 v 3 - Y v Yi~ 1 Y v 



T v 2 
Vl 



(7) 



Since in such "1-2-3" seesaw lepton number is ungauged, there is a physical Goldstone boson resulting 
from its spontaneous breakdown, namely the Majoron [26 27]. It is often argued that, due to quantum 
gravity effects the associated Majoron will pick up a mass. It has been shown that, a keV range Majoron 



can provide the observed dark matter of the Universe 1 28 1 and be detected through its X-ray gamma line 
searches [29|. If B-L is gauged [ 30 1 the Majoron is absorbed as the longitudinal mode of a new neutral 
gauge boson. 



2.2 Low-scale seesaw mechanisms 

A distinguishing feature of the seesaw mechanism as proposed in Ref. [23 25 26 1 and other presenta- 
tions p4) is that it is formulated in terms of the standard SU(3) S3 SU(2) £g> U(l) SM gauge group. 
The higher generality implies that the number of "right-handed" neutrinos is totally arbitrary since, being 
gauge singlets, they carry no anomaly. New important features may emerge when the seesaw is realized 
with non-minimal lepton content, opening the door to the possibility of low-scale seesaw mechanisms, 
such as the inverse seesaw (3TJ. 



2.2.1 Inverse seesaw mechanisms 

The model adds a pair of two-component SU(3) ® SU(2) (g> U(l) singlet leptons, uf, Si, to each SM 
generation i running over 1. 2, 3. In the basis v, v c : , S, the neutral leptons mass matrix M. v is 9 x 9, i.e. 




M v = Y v v 2 M l \, (8) 



with /i <C Y v v% <C M, where Y v and M are arbitrary 3x3 complex Yukawa matrices, while // is complex 
symmetric, due to the Pauli principle. In such a scheme the three light neutrino masses are determined 
from the effective 3x3 neutrino mass matrix 

m v pb v 2 2 Yi?M T ~ 1 iiM- 1 Y v . (9) 
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Fig. 3 Inverse seesaw mechanism. 



The mass generation is illustrated in Fig. [3] Notice that as/i->0 all neutrinos become massless and lepton 
number symmetry is restored. The entry fi may be proportional to the vev of an SU (3) ® SU{2) ® £7(1) 
singlet scalar, in which case the model contains a singlet Majoron [32] which may provide an invisible 
Higgs boson decay channel [33]. In such schemes one must take into account the existence of sizeable 
invisible Higgs boson decay channels in the analysis of experimental data on Higgs searches J34][35] | . 



2.2.2 Linear seesaw mechanism 

We now turn to a low-energy seesaw mechanism with gauged B-L, originally suggested in the framework 
of dynamical left-right symmetry [36, 37 1 and more recently in an SO{10) model with broken D-parity in 
which only an Abelian factor survives at low energies J38J. In addition to the three left- and right-handed 
neutrinos the model contains three sequential gauge singlets So, with the following mass matrix 

Y v (*) 

M v = I Yj » F( XR ) I (10) 



F t (xl) F t ( X 




in the basis Vl, v c l , Sl- The zeros along the diagonal, in the v^-vl and v c L -v c L entries, are due to the fact 
that there is no 126. The resulting neutrino mass is 



Y V {FF- X f + (FF- 1 )Y, 1 



(11) 



where M un jf is the unification scale, F and F denote independent Yukawa coupling combinations of the 
SiL- One can see that the neutrino mass is suppressed by the unification scale M nn a instead of the B-L 
breaking scale. Note that, in contrast to other seesaw schemes, this one is linear in the Dirac Yukawa 
couplings Y v , as illustrated in Fig. [4] It is rather remarkable that one can indeed take the B-L scale as 



$ \ /' XR 



X X 



I I I 
\XL J XR J $ 
I I I 
-* 4 * ft * 



3 



Fig. 4 Linear seesaw mechanism. 



low as TeV without generating inconsistencies with gauge coupling unification [38 1. The light neutral 
gauge boson may be searched directly at the Large Hadron Collider (LHC) or through precision studies of 
low-energy neutrino-electron scattering |39). 
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2.3 Radiative models of neutrino mass 

In addition to the above low-scale seesaw schemes there is a variety of other models of neutrino mass 
where the operator O a b is induced from physics at accessible scales, TeV or less. The first possibility is 
that neutrino masses are induced by calculable radiative corrections pOpT) , for instance, as illustrated in 
Fig. [5] Up to a logarithmic factor one has, schematically, 



V167T 2 / (m k y 

in the limit where the doubly-charged scalar k is much heavier than the singly charged one. Here I denotes 
a charged lepton, / and h are their Yukawa coupling matrices and YJ denotes the SM Higgs Yukawa 
couplings to charged leptons and (a) is an SU(3) (8 SU(2) <g> £7(1) singlet vev introduced in Ref. |42|. 
The smallness of the neutrino mass arises from the presence of a product of five small Yukawas and the 
appearance of the two-loop factor. A special feature of the model is that, thanks to the anti-symmetry of 
the / Yukawa coupling matrix, one of the neutrinos is massless. 



Yl r R £ i R i L v« c 

Fig. 5 Two-loop origin for neutrino mass. 



2.4 Supersymmetric neutrino masses 

An interesting alternative are models where low energy supersymmetry is the origin of neutrino mass p3) 
through the breaking of the so-called R parity. This could arise spontaneously, driven by a nonzero vev 
of an SU(3) ® SU(2) ® U(l) singlet sneutrino |44|]46| . This way we are led to the minimal way to 
include neutrino masses into the MSSM, which we take as reference model, with effective bilinear R parity 
violation |47|. The neutrino mass generation scenario is hybrid, with one scale generated at tree level by 
the mixing of neutralinos and neutrinos, and the other induced by "calculable" loop corrections J48]|49) . 
The neutrino mass spectrum naturally follows a normal hierarchy, with the atmospheric scale generated at 




Fig. 6 Loop origin of solar mass scale. Atmospheric scale arises from tree-level neutralino exchange, 
the tree level with the solar mass scale arising from calculable loops, as indicated in Fig. [6] 

3 Prototype flavour model with tetrahedral symmetry 

We now turn to models incorporating flavour symmetries, starting with the BMV (Babu-Ma-Valle) model p3] . 
The usual quark, lepton, and Higgs superfields transform under A4 as follows: 
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The following heavy SU(2) singlet quark, lepton, and Higgs superfields are also added 





u 


U c 
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with a; 3 = 1 and 1 + ui + uj 2 = 0. The superpotential of this model is then given by 

W = MuUiU,, 



fuQiUtfa + K^Uiu'xk + M D D t D c t + fdQiD'fa + h% k Dd c ,Xk 

1 



+ M E EiE? + f e L t E;4i + h% k E t e c 3 Xk + \M N NfN? + f N UN^ 2 + ^ifo 
+ \M x XiXi + h x XiX2X3, 



(13) 



with the usual assignment of R parity to distinguish between the Higgs superfields, i.e. an d Xi> from 
the quark and lepton superfields. The terms XiNjNg, etc. are forbidden by the Z%. However, Z 3 can break 
explicitly but softly, by M x ^ 0. The scalar potential involving Xi is given by 



V = \M xX i + h xX 2X3\ 2 + \M xX 2 + h xX3 Xi\ 2 + \M xX3 + KXiXi? , 
which has the supersymmetric solution (V = 0) 

(Xi) = (X2) = (X3)=u = -M x /h x , 



(14) 



(15) 



so that the breaking of A 4 at the high scale M x does not break the supersymmetry. Consider now the 6 x 6 
Dirac mass matrix linking (e i7 E{) to (e|, Ej). 



M 



eE 












feVl 




















feVl 




















feVl 


hi u 


h\u 


h\u 


M E 








h\u 




h\uj 2 u 





M E 





h\u 


h%uj 2 u 


h\iou 








M E 



(16) 



where vi = (cj)f) with similar forms for the quark mass matrices. The reduced 3x3 charged leptons mass 
matrix is then 



M e =U L 



hi' 






h e 2 ' 






hf 



V3feVlU 



Mi 



(17) 



where hf = hf[l + (hfu) 2 /Mf]" 1 / 2 and 



Ur 



1 



l l l 

1 UJ 
,2 



2 



(18) 



This shows how charged-lepton masses are allowed to be all different, despite the presence of the A 4 
symmetry, because there are three inequivalent one-dimensional representations. Clearly, the up and down 
quark mass matrices are obtained in the same way, both are diagonalized by Ul, so that the charged-current 
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mixing matrix Vckm is the identity matrix. CKM angles may be generated from corrections associated to 
the structure of the soft supersymmetry breaking sector |50 51 1. The 6x6 Majorana neutrino mass matrix 
is given by 



M v 



N 



U L f N v 2 
Ulf N v 2 M N 



(19) 



in the basis (v ei v T ,N£, 7V|, N§) and v 2 = {(j^)- The effective (v e , v T ) mass matrix becomes 



M 



N 



M 



N 



1 
1 
1 



(20) 



showing that neutrino masses are degenerate at this stage. Consider now the above as coming from an 
effective dimension-five operator 



M N 



(21) 



where A ee = A MT = A rM = 1 and all other A's are zero, at some high scale. As we come down to the 



electroweak scale, Eq. (14) is corrected by the wave-function renormalizations of v e 



and v r , as well 



as the corresponding vertex renormalizations, lifting the neutrino degeneracy due to the different charged- 
lepton masses. In order to obtain a pattern suitable for explaining current neutrino oscillation data we 
assume the presence of radiative corrections associated to a general slepton mass matrix in softly broken 
supersymmetry (521. Given the structure of Ay at the high scale, its low scale form is fixed to first order as 



1 + 2S ee 
5 ef i + 8 eT 
8 eli + 8 eT 1 



Sen + 



f-5 



8^ + S T 



SeT 

'jlljU "T" 8 rT 
28^ T 



(22) 



where we have assumed all parameters to be real as a first approximation. [The above matrix is obtained by 
multiplying that of Eq. (13) on the left and on the right by all possible v. L — > Vj transitions.] Let us rewrite 



the above with S = 8^ + S TT — 2S^ T , S = 25 /iT , 8 
1 



- W 2 - S tt/2 - 8„ T , and 8" 



A, 



■S + 25 + 28' 
8" 
5" 



1 



5" 8" 
5 1 + S a + 8 
8 + 8 8 



so that the exact eigenvectors and eigenvalues are easily obtained: 



and 















Ai - 


H 
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A 3 = 


-1 



cos 6 sin 9/V2 sin0/\/2 
-sin cos^/a/2 cos0/-\/2 
-1/V2 1/^2 



v.. 



leading to 



(23) 



(24) 



(25) 
(26) 
(27) 



sin 2 2Q atm = 1, tan 2 1 



'sol = 



5 n2 + s ,2 _ 5 i^ 5 a + 2< 5" 2 ' 



(28) 
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if 8' < and \6"/6'\ = 1.7. Assuming 8" 2 /8' 2 and 8', 8" « S, we now have 



Am 2 32 



ASml, Am?, ~ 



28" 2 m 2 



(29) 



where mo is the common mass of all 3 neutrinos. This provides a satisfactory first-order description of 
present neutrino-oscillation data. With f7 e 3 — there is no CP violation in neutrino oscillations. However, 
if we assume complex A,j then one has one CP phase which cannot be rotated away. Without loss of 
generality, we now rewrite \j as 



K — 



1 + 25 + 28' 8" 5"* 
8" 8 1 + 8 

8"* 1 + 8 8 



(30) 



where we have redefined 1 + So as 1, and 8, 8' are real. Assuming that 8', ReS" and (Im8") 2 /S are all 
much smaller than <5, one can diagonalize this mass matrix approximately, 



v , (ImS") 2 
S+ ^28—> 3 



ReS". 



(31) 



obtaining that U e s is imaginary and CP violation in neutrino oscillations is predicted to be maximal. There 
is also an interesting relationship, i. e. 



'12 



-1 2 



\u, 



ReS" 



(32) 



indicating that \U e ^\ is naturally of the order | Ara^/Aroy 1 / 2 

lU^^OdAmlJAm 2 ^/ 2 ) 



in the limit 8', 8" < 8. 

In FigjTjwe plot the maximum achievable value of Am 2 tla against the overall neutrino mass scale mo- 
The value of mo is subject to an upper bound given by (/3/3)oi> experiments ]53|54[ and cosmology J55p6) . 
In order to get large enough neutrino mass splittings to account for current oscillation data we need 



m > 0.4 eV 



(33) 



leaving a relatively small room for the value of mo. Note also that the mass splittings are related to the 
parameters 8, 8' and 8", and these are increasing functions of the slepton mixings and also mass splittings. 
This is potentially in conflict with the restrictions from lepton flavour violation searches which push the 
spectrum toward mass degeneracy and small mixings. However one can show that viable spectra do exist. 



4 Quarks, non-abelian discrete flavour symmetries and unification 

In the original BMV model fl5| as well as in the subsequent Altarelli-Feruglio model JT7| the CKM mixing 
matrix was predicted to be the identity, which provides indeed a good first order approximation. However, 
realistic quark mixings require either renormalization effects 1 50|51 1, or suitable model extensions, e.g (57) 
and ||58). 

The largest angle in the CKM matrix is the Cabibbo angle governing the mixing between first and 
second generations, about Ac ~ 0.22. Mixing angles between the first/third and second/third families are 
about A c and A c , respectively. This suggests that first and second quark families belong to a doublet of a 
family symmetry, instead of singlets or triplets. 

A non trivial extension of the Altarelli-Feruglio model for quarks was given in J59) where A4 is extended 
to its double covering T". The main advantage of this group is that it contains doublet irreps besides the 
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Fig. 7 Maximum atmospheric mass squared difference versus mo (light shaded histogram). The horizontal band is 
the 3a allowed region for AmJ tm by current MINOS/T2K data. 



triplet. This feature of T' is suitable for the quarks since the first two generations can be assigned to doublet 
irreps, while the third generation belongs to a singlet. 

Another group which is interesting in order to include quarks in a TBM pattern for leptons is the permu- 
tation group of four objects 54. Lam |60| noted that the mininmal group for TBM is S4 and Refs. ]6T][62) 
have studied a model based on S4. In Ref. J63) it has been shown that the Cabibbo angle can be predicted 
using the dihedral flavour symmetry D n . 

However in general it seems that there is not yet a simple and elegant framework that can explain at the 
same time neutrinos and quarks. Most of the models considered for quarks contain many flavon fields and 
extra ad hoc abelian symmetries. 

Recently there has been a lot of effort in studying the possibility to embed the TBM ansatz within a 
grand unified (GUT) framework. The most popular unifying groups considered were SU(5), Pati-Salam 
SU C (4:) x 5C/l(2) x SUr(2) and 50(10). These can be separated into two classes with respect to TBM. 
In the SU(5) framework it is easier to obtain the TBM pattern than in the case of Pati-Salam or 50(10), 
since right-handed neutrinos transform trivially under the gauge group. For some example of extensions 



of SU(5) with discrete flavour symmetries see for instance Ref. |64-70|, with Pati-Salam ||7lj|72| and 



50(10) [73-82]. 



5C/(5) models 

As an illustrative example here we consider the model studied in |64j. In typical SU (5) GUT models, 
the 5 contains the lepton doublet L and the three (colored) right-handed down type quarks d c for each 
family. Whereas the 10 contains the right handed charged leptons, the three right handed up type quarks 
u c and the three quark doublets Q for each family. Typically in ^-based models the three lepton doublets 
are assigned as triplets of A4, while the right-handed charged leptons are assigned to A4 singlets. This 
means that we should choose 

5-3, 10~1,1',1" (34) 
which implies the following quark assignments 

d?~3, (Ui,di)~ 1, 1', 1", (35) 

One assumes three Higgs doublets 5# — 3 in the down/lepton sector 

hi \ / Vi 

U, = i ha \ U L • v 2 I (36) 
h 3 I \ v 3 
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where hi are Yukawa couplings and Vi = (5 H . ) |j For the up-quark sector one assumes only two Higgs 
doublets, namely 5h ~ 1', 1", so that the up quark mass matrix has the form 

/ ft ft \ 

M u = \ H2 m 2 (37) 
\ /i 3 m 3 J 

where 7712 and ^3 arise from the vev of 5h ~ l' and 7713 and /12 come from the vev of 5h ~ l". 

In Ref. J64) it has been shown how the CKM mixings V^ s , K( S and V u t, as well as the CP phase can 
be obtained in this model. In total the model has 10 free parameters in the quark sector for 10 observable 
(6 masses, three mixing angles and one phase). Therefore it is clear that in principle such a model can 
fit realistic quark masses and mixings. Indeed this has been shown, however no light can shed into the 
structure of quark masses and mixings. Moreover we has the usual SU(5) relations m T — nib, "V = Tn a 
and m e = at the GUT scale. While the first relation is in good agreement with data, the last two are not. 
To decouple the charged lepton sector from the down quark sector the usual strategy is to use bigger ST/ (5) 
representations in the scalar sector, like the 45. In the simplest scenario the gauge coupling unification has 
not been considered. 



50(10) models 

The situation in 50(10) with type-I seesaw is much more complicated than in 5t/(5). Indeed the main 
problem is that in 5O(10) TBM requires to distinguish the Dirac neutrino Yukawa coupling from that of 
the up-quark. In particular the former must be either proportional to the identity or given by eq.([T]l in the 
basis where charged leptons are diagonal. In contrast the up quarks must be strongly hierarchical, namely 



M u ~ x x x , m D ~ 1 , (38) 



X 


X 


X 


X 


X 


X 


X 


X 


1 




where x indicate small entries much smaller then one. We consider for simplicity the case where the 
Dirac neutrino mass matrix is proportional to the identity. As already mentioned, in 5O(10) all matter 
fields (including right-handed neutrinos) belong to only one multiplet, the spinorial 16. In renormalizable 
50(10) models only three types of Yukawa contractions are possible, namely 16-16 = 10+120+126. Two 
of them Y"i u and Yj^ are symmetric and Y 12 o is antisymmetric. It is well known that their contributions to 
up quark and Dirac neutrino mass matrices are given as 

M u = M v = Y 10 i; 5 , M" = Yj^r 5 , M v = -3Y^r 5 , M u = Y 12Q u 45 , M v = Y 120 u 5 (39) 

where v a , u a ,r a are the vevs of 10, 120 and 126 respectively, and the lower indexes a = 5, 45 are SU(5) 
components. Note that the 10 gives equal contributions to the Dirac and up-quark mass matrices, while the 
126 gives contributions to the Dirac and up-quark mass matrices that are proportional to each other. There- 
fore if we take hierarchical Yxo an d ^126 trie Dirac neutrino mass matrix will be hierarchical, in contrast 
with TBM requirement given in eq. p8j ). The 120 gives different contributions to the up quarks and Dirac 
neutrino mass matrices, proportional to u 45 and u 5 , respectively. Therefore one can in principle distinguish 
between Dirac neutrino and up quarks by means of the 120. However the coupling Y120 is antisymmetric 
giving two degenerate eigenvalues and zero determinant. Hence it is not possible to obtain a hierarchical 



Yukawa coupling with 120. We conclude that it is not possible to obtain the TBM mixing pattern eq.(38 1 
within a simple renormalizable 5O(10) framework with type-I seesaw mechanism as described above. 
This problem has been circumvented in two ways: 

• assuming type-II dominant with respect to type-I seesaw; 



^ Since in minimal SU(5) we have that M; = Mj it is clear that one can have TBM mixing in the limit v\ = V2 = v 3 . 
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introducing non-renormalizable operators. 



The idea of type-II dominance was suggested in Ref. |83| in the context of an 5*0(10) model with 
quasi-degenerate neutrinos. A similar scenario has been used in Ref. |79) to accommodate the TBM 
mixing pattern. The idea is that in 50(10) the type-II seesaw arises from the coupling with 126 scalar and 
is proportional to its 15 component under SU (5). It is well known that the 15 does not give contributions 
to the Dirac fermion masses, see eqs.([39]>. Therefore assuming type-II dominance neutrino and charged 
fermion Yukawa couplings become unrelated and we can easily obtain TBM mixing pattern. 

In fact we can take the Yukawa coupling Yi 2 q with TBM form given in eq.([T]i, while Dirac neutrino 
mass matrix and up quark mass matrix can be taken hierarchical. This is no longer a problem since type-I 
seesaw contribution is assumed to be negligible with respect to type-II contribution, giving only deviations 
from the TBM pattern that can generate a sizable #13 angle. 

The second possibility to address the problem of having a TBM mixing patter within the 50(10) type-I 
seesaw mechanism is to use non-renormalizable operators. For instance one can use the dimension five 
operator^] 

/i 16 16 120// 45//. (40) 



This gives a contribution to the up-quark mass matrix but not to the Dirac neutrino 1 75 84 85 1 . This can be 



seen in more details as follows. The 45// can take vev in its singlet ls<y(5) component called X-directiorl 5 
or along the adjoint 24.su (5) component, that is the hypercharge y-direction (see for instance |86)). We 
indicate their vevs as 

h = (W(5)>, hi = (24^(5)). (41) 

The SU(5) components of the 120// of 5O(10) that contain SU(2) doublet (giving rise to the Dirac 

masses terms for the fermions) are the 455/7(5), 45g//(5), §su(5) an d 5su(5) representations. Denoting 
their vevs as 

a 5 = (5g//(5)), a§ = (5,5(7(5)), a 45 = (45 S£ /( 5 )), a 4 - 5 = (45 S£ /( 5 )), (42) 
we find that 

M u = h a 45 (h -4b 2i )-h T a i5 (b 1 + b 2A ), (43) 

M v = 5ha 5 bi - h T a 5 (-36i - 3b 24 ), (44) 

M d = h{a^ + a IS )(-3b 1 +2b 2 4)- h T {a^ + a TS )(b 1 +b 2 4), (45) 

Mj = h{a z - 3fl45)(-36i - 36 24 ) - /i T (% - Za^h + 6b 24 ). (46) 

where h is an arbitrary matrix. Note that if we set 624 = 0, M u is proportional to h — h T which is 
antisymmetric. In contrast, if we set bi — and b 2 4 7^ then M u oc Ah + h T which is an arbitrary matrix. 
When 05 = the operator 16 16 120// 45// contributes to Y u , Yd, Y e but not to Y v , so we can distinguish 
between Dirac neutrino and up-quark mass matrices. 



Recent observations of a relatively large reactor angle [6-8] pose the questions as to whether the TBM 
pattern is indeed a good starting point to describe neutrino flavour mixing. The difficulty encountered in 
5O(10) based GUT models suggests us to discard the TBM ansatz as starting point. However discrete non- 
Abelian flavour symmetries appear to be better candidates for a family symmetry. Here to give an explicit 
example of such possibility, briefly presenting a model given in Ref. [87| based on 5O(10) x D 3 | 
The group D3 (isomorphic to 53-the group of permutations of three objects) contains three irreducible 



This operator can be obtained by integrating out a couple 16 x — 16 x with renormalizable couplings 16 16 x 120// > 16 16 x 45// . 
5 This is the extra 1/(1) contained in 5O(10) D 51/(5) X J7x(l). 
" An extra U(X) family symmetry is also required. 
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representations, one symmetric singlet, one antisymmetric singlet and one doublet. The third generation is 
assigned to the antisymmetric singlet, I63 ~ 1a while the first two families are assigned to a doublet of 
D3. The only renormalizable coupling is for the third generation, with first and second generation masses 
arising only from next to leading order contributions. Extra messenger fields are introduced in order to 
make the Lagrangian renormalizable. All charged fermion Yukawa couplings Y u , Y^, Y( and Y v have 
Fritzsch texture [88|, namely 




Y f oc x x x . (47) 



Such a model is that it contains only seven real free parameters plus three complex phases for thirteen 
observable charged fermion masses and mixing angles. While in the neutrino sector there are other four 
free parameters, giving in total 14 free parameters. Hence this model is highly predictive. 

5 Dark matter and flavour symmetry 

Deciphering the nature of Dark Matter (DM) constitutes one of the main challenges in cosmology since 
decades. Some recent direct and indirect DM detection experiments have given tantalizing hints in favour 
of a light WIMP-like DM particle |89}j93) feeding the hopes of an imminent detection. 

An interesting idea investigated recently has been to link neutrino mass generation to dark matter, two 
seemingly unrelated problems, into a single framework. This is not only theoretically appealing, but also 
may bring us new insights on both issues. 

Among the requirements a viable DM candidate must obey, stability has traditionally been ensured by 
imposing of a stabilizing parity in an ad hoc way. It would be clearly appealing to obtain stability in a 
theoretically natural way. This has motivated attempts such as gauged U(1)b-l Q94L gauged discrete 
symmetries ]95) as well as the recently proposed discrete dark matter mechanism (DDM) [96-99], where 
stability arises as a remnant of a suitable flavour symmetry]^] 

The interplay between decaying dark matter and non-Abelian discrete flavour symmetries has been 
considered in a number of papers; for instance, in | |103f 105 1 non-Abelian discrete symmetries prohibit 
operators that may induce too fast dark matter decay; in | |106| a non-Abelian discrete symmetry (not a 
flavour symmetry) has been used to stabilize the scalar DM candidate (similar to what has been discussed 
in the inert scalar models p07| ). However, these models differ substantially from what was proposed 
in (55). 

Here we describe the original suggestion. Consider the group of the even permutations of four objects 
A4. It has one triplet and three singlet irreducible representations, denoted 3 and 1, 1', 1" respectively. 
A4 can be broken spontaneously to one of its Z2 subgroups. Two of the components of any A4 triplet are 
odd under such a parity, while the A4 singlet representation is even. This residual Z2 parity can be used to 
stabilize the DM which, in this case, must belong to an A 4 triplet representation, taken as an SU (2)^ scalar 
Higgs doublet, rj ~ 3 [96-99]. Assuming that the lepton doublets Li are singlets of A4 while right-handed 
neutrinos transform as A4 triplets N ~ 3, the contraction rules imply that 77 couples only to Higgses and 
heavy right-handed neutrinos L{ Nfj. In this case 77 and N have even as well as odd-components while Li 
are even so that Li N 77 interaction preserves the Z2 parity. Invariance under Z2 implies that N components 
odd under Z 2 are not mixed with the Z 2 -even light neutrinos z^. This forbids the decay of the lightest Z 2 - 
odd component of 77 to light neutrinos through the heavy right handed neutrinos, ensuring DM stability. In 
this framework DM has quartic couplings with the SM Higgs doublet as is the Higgs portal DM scenario, 
and has been shown to have the correct relic density |108) , with annihilation and co-annihilation of DM 
into SM particles (fermions and bosons), see Fig.|8]for the results. Note that assigning the three left-handed 
leptons to a flavour-triplet implies that the "would-be" DM candidate decays very fast into light leptons, 
through the contraction of the triplet representations, see general discussion in ref. [ 1 00 1 . This problem 



For other flavour models with DM candidates see |97|l00]{l04) 
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M DM [GeVJ M DM [GeVJ 

Fig. 8 Left plot: Regions in (Mdm - Mb) plane (DM mass-lightest Higgs boson mass) allowed by collider con- 
straints and leading to an adequate DM relic abundance. Right plot: Spin-independent DM scattering cross section 
off-protons as a function of the dark matter mass. From Ref. |98|. 



has been considered by Eby and Framptom [ 109 1 using a T' flavour symmetry. While the suggested model 
has the merit of incorporating quarks non-trivially, it requires an "external" Z2 asymmetry in order to 
stabilize dark matter. In fact this observation lead ref. ]1 10[ to claim that a successful realization of the 
DDM scenario requires the lepton doublets to be in three inequivalent singlet representations of the flavour 
group. Recently Ref. [ 108 1 has given an explicit example of a model based on a A (54) flavour symmetry 

0.20 1 1 „ _ mi . , 




'0.20 0.25 0.30 0.35 0.40 

sln 2 Si2 

Fig. 9 The shaded (yellow) curved band gives the predicted correlation between solar and reactor angles at two-sigma 
for normal hierarchy. The solid (black) lines give the global best fit from Ref. 1 12 J. The dashed lines correspond to the 
central values of the recent reactor measurements |6- 8|, along with the corresponding two-sigma bands. 

in which left-handed leptons are assigned to non-trivial representations of the flavour group, with a viable 
stable dark matter particle and a nontrivial inclusion of quarks. In contrast to the simplest "flavour-blind" 
inert dark matter scheme JTTTJ, such a model implies non-trivial restrictions and/or correlations amongst 
the neutrino oscillation parameters, consistent with the recent reactor angle, see Fig. [9] Although neutrino 
mixing parameters in the lepton mixing matrix are not strictly predicted, as seen in Fig. there is a 
non-trivial correlation between the reactor and the atmospheric angle. While the solar angle is clearly 
unconstrained and can take all the values within in the experimental limits, a nontrivial correlation exists 
with the reactor mixing angle, indicated by the bands in Fig. [9] These correspond to two and 3<r regions 
corresponding to the global oscillation fit in Ref. fT2| . The horizontal lines give the best global fit value 
and the recent best fit values obtained in Daya-Bay and RENO reactors j7j[8]. 



6 Neutrinoless double beta decay 

This is the process par excellence which allows us ways to test the fundamental nature - Dirac or Majorana 
- of neutrinos in a model-independent way, i. e. irrespective of which mechanism generates neutrino 
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masses and irrespective of which mechanism induces Of/3/3 (neutrinoless double beta decay) The basic 
connection is given by the black box theorem | 1 16 1 17 1 illustrated in Fig. 10 Moreover, Of /3/3 receives a 



w 




w 



Fig. 10 The observation of Of/3/3 implies the Majorana nature of a neutrino |116| . 

contribution from the tree-level exchange of Majorana neutrinos, whose amplitude, illustrated in Fig.fTTj 
is proportional to an effective mass parameter which, in contrast to neutrino oscillations, is sensitive also 
to the absolute scale of neutrino masses, which is independently tested also in searches for tritium beta 
decay (TT8][n9j and cosmology j55]|56). 

In addition, this amplitude can bring complementary information on the underlying flavour structure as 
revealed, say, in neutrino oscillation searches^] 

As we saw in Sec.|3|the BMV model [ 15] implies a lower bound on the absolute neutrino mass m u > 0.4 
eV and therefore wilibe tested fairly soon in Of /3/3 searches. 

On the other hand, even if neutrinos are non-degenerate, many of the models based on non-Abelian dis- 
crete flavour symmetries are characterized by a specific (complex) relation between neutrino mass eigen- 
values, leading to mass sum rules (MSR). Such MSR lead in most of the cases to lower bounds for the 
neutrinoless double beta amplitude parameter, depending of the type of the specific MSR. The following 
types of mass relations hold: 

(48) 
(49) 
(50) 
(51) 



Here m\ = m\ denote neutrino mass eigenvalues, up to a Majorana phase factor, while \ ar, d £ are & ee 
parameters which specify the model, taken to be positive without loss of generality. 

We first consider the amplitude for neutrinoless double-/3 decay within a flavour-generic scheme. The 
effective neutrino mass parameter \m ee \ determining the Of/3/3 decay amplitude, as a function of the light- 
est neutrino mass is given in Fig. [TT] As is well-known, by varying the neutrino oscillation parameters 
in their allowed ranges Jl2} one obtains two types of relatively broad bands in the (|m ee |, m^„ ht ) plane 
il and 




corresponding to normal and inverse hierarchy spectra, as shown in Fig. 1 1 For such "flavour generic" 
case one finds a lower bound for the neutrinoless double-/3 decay effective mass parameter \m ee \ only in 
the case of inverse mass hierarchy. Indeed, thanks to the possibility of destructive interference among the 



light neutrinos no lower bound holds for the case of normal neutrino mass hierarchy 1 1 1 2 121 122 1. 

Let us now turn to the case where "flavoured" case where MSR relations like (A),(B),(C) and (D) hold. 
As mentioned above these can be obtained in flavour models where the neutrino mass matrix depends 
only on two independent free parameters, so that the resulting mixing angles are fixed, as in the case of 
tri-bimaximal or bimaximal mixing patterns. 



CP and electromagnetic properties of neutrinos |1 12||1 15| are also sensitive to the fundamental nature of neutrinos, though 
experimental prospects are far less clear than those for Of/: 
9 



Subject, of course, to nuclear matrix element uncertainties 120 . 
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Fig. 11 Allowed range of (|m ee |) as a function of the lightest neutrino mass for the TBM mixing pattern (red and 
green bands for NH and IH respectively) and for the full 3a C.L. ranges of oscillation parameters from 1 12 1 (gray and 
blue bands for NH and IH respectively). 



For definiteness here we focus on the case of tri-bimaximal neutrino mixing pattern. Taking into account 
that corrections from higher dimensional operators and/or from the charged lepton sector can yield (9 13 ^ 0, 
here we retain the TBM approximation as a useful starting point to obtain our MSR relations. However, 
when evaluating a lower bound on the effective neutrino mass parameter |m ee | determining the Qvf3j3 decay 
amplitude, we include explicitly the effects of non-vanishing #13, by taking the 3 a oscillation parameter 
values determined in Ref. fTZ] . 

This way one obtains lower limits of \m ee | corresponding to different integer choices of (x, £) between 1 



and 3 and for each of the four MSR considered in eqs. (48 i-(51 1, for both normal and inverted hierarchies. 
These results are summarized in Tab. I of Ref. |123) . A large class of non-Abelian flavour symmetry 
models discussed in in the literature is covered, see for example, references [17-19,57,58,61 62,65-67, 
|75l [T24HT38l[T38 ri T48| - 

Some comments are in order. First let us consider the effect of a possible non-zero effect of #13 as 
indicated by recent experiments [5 6] as well as global neutrino oscillation fits |12 149) . In Fig 12 we 
show the prediction for |m ee | as function of mught obtained from the MSR 3^/rfi2 + 3y / m3 = ^Jm\ (right 
panel) and IJrn^ + Jm 3 = J mi (left panel). For the red bands we assumed the TBM values of the 



oscillation parameters (implying (9 13 = 0) while the yellow bands corresponds to the same MSR, but now 
varying the values of #13, 6*23 and #12 within their allowed 3a C.L. interval. By looking at the left panel in 
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Fig. 12 ( | m ee | ) as a function of the lightest neutrino mass corresponding to the mass sum-rule 2^/7712 + s/mi = %/mi 
|150| (left) and 3^/7712 + 3-^/772,3 = *Jm\ (right). The red bands correspond to the TBM mixing pattern, while the 
yellow bands correspond to the same MSR, but now varying the values of #13, 6*23 and #12 to 3a C.L. range. 



Fig. 12 one sees that, indeed, the Qvf3(3 lower bound is sensitive to the value of #13. 
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One also finds that, as expected on general grounds, all inverse hierarchy schemes corresponding to 
various choices of (x, £) within sum-rules A-D have a lower bound for the parameter |m ee |. However, the 
numerical value obtained depends on the MSR scheme, signaling that not all values within the correspond- 
ing band in Fig.[TT]are covered for a given flavour symmetry structure. 



7 Lepton flavour violation and flavour symmetry 

Flavour violation is required to account for the current neutrino oscillation data. It should, however, make 
its appearance also in other sectors, inducing rare proce sses invol ving charged leptons, whose strength 
is not suppressed by the smallness of neutrino masses 1 151 152| [*°| For example, in the presence of 
supersymmetry, the lepton flavour violation required to account for oscillation data in high-scale seesaw 
schemes induces decays such as fi~ — > e~7 and flavour violating tau decays (Fig. 13 1 as well as nuclear 
fi~ — e~ conversion (Fig. 14 1 as a result of the exchange of supersymmetric leptons. The existence of 
such loop effects, illustrated in Fig. 13 has been known for a while |51 151 1. These lepton flavour 





Fig. 13 Feynman diagrams for l i — > I - 7 involving chargino/neutralinos and sneutrino/charged slepton exchange. 



r 



7 \ q 

"TJucleon" 



c_ _Nucleon_ ^ 
" Nucleus 



(a) (b) 

Fig. 14 Contributions to the nuclear /i~ — e~ conversion: (a) long-distance and (b) short-distance. 



violation processes can have interesting rates, which depend not only on the seesaw mechanism, but also 
on the details of supersymmetry breaking and on a possible theory of flavour. The resulting lepton flavour 



violation rates will be accessible to the upcoming generation of experiments [ 155 - 157 1. 



7.1 Non-supersymmetric low-scale seesaw models and flavour 

Low energy seesaw schemes such as the inverse or linear seesaw mechanisms generate neutrino masses 
from fermion messengers (right-handed neutrinos) at the TeV scale pT|[36-38|. These are potentially 
accessible to the LHC, especially in the presence of a new gauge boson "portal" associated, for example, 
to left-right symmetry |158||159) . 

Within low-scale seesaw mechanisms lepton flavour violation and/or CP violating effects arise at the 
one-loop level from the exchange of relatively light neutral heavy leptons. Their strength is not sup- 
pressed by the smallness of neutrino masses P51||152| so the resulting lepton flavour violation effects 



10 



Similar results hold also for leptonic CP violation 153|l54[ 
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are potentially large even in the absence of supersymmetry 1 152 - 154) 1 160 161| and/or extended gauge 
structure fl58l[l59l 



In the inverse and linear seesaw models proposed in [ 162 1, the neutrino mass matrix is a 9 x 9 symmetric 
matrix. It is diagonalized by a unitary matrix U a p, a 1 (3 = 1...9, leading to three light Majorana eigenstates 
Vi with i = l,2,3 and six heavy ones Nj with j = 4, .., 9. The effective charged current weak interaction 
is characterized by a rectangular lepton mixing matrix Ki a [25 1, 



C CC = ~K ia Li^(l + l5 )N a W, 



(52) 



where i = 1, 2, 3 denote the left-handed charged leptons and a the neutrals. The contribution to the decay 
h —> Ijl arises at one loop from the exchanges of the six heavy right-handed Majorana neutrinos Nj which 
couple sub-dominantly to the charged leptons. The well-known one-loop contribution to this branching 
ratio is given by ]161[ 



Br(k -> l n ) 



1 



a 3 sw mf. 
256tt 2 NT*, T h 



\Gij\ 2 



(53) 



where 



Gij —J2k-- 

GJx) 



2x 3 +5x 2 -x 



Ml 
3x 3 



(54) 



i(l-x 3 ) 



2(l-z) 



In a; 



We note that, thanks to the admixture of the TeV neutral leptons in the charged current weak interaction, 
this branching ratio can be sizeable [152|. Similar results hold for a class of LFV processes, including 
nuclear mu-e conversion [ 163 1 whose expected rates are strongly correlated to those of // — » &y, see Fig. [18] 
As an example we now consider the low-scale seesaw TBM model given in |162| . The simplicity of their 
mass matrices, which are expressed in terms of very few parameters, makes such a models especially 
restrictive and this has an impact in the expected pattern of LFV decays. In contrast to the general case 
considered in 1 163 164|, in 1 162 1 one can easily display the dependence of the fi — > branching ratio on 
the new physics scale represented by the parameters M ~ TeV and the parameters /i or vl characterizing 
the low-scale violation of lepton number. This is illustrated in Fig. [T5| 




Fig. 15 Br(/j, — > erf) versus the lepton number violation scale: n for the inverse seesaw (red color), and Vl for the 
linear seesaw (blue color). Here M is fixed as M = 100 GeV (continuous line), M — 200 GeV (dashed line) and 
M = 1000 GeV (dot-dashed line). 



1 1 In type-I seesaw schemes the processes I, 
not enough to give large rates since the messenger scale Mr characterizing lepton number violation is too high 



lj 7 are enhanced due to a breakdown of the GIM mechanism. However this is 
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Note also that, in contrast to a generic inverse or linear seesaw model, such a A,4-based models the 
structure of the matrix Gij is completely fixed, and this leads to predictions for ratios of lepton flavour 
violation branching ratios such as 

Bt(t — > ey) / m r \ 5 r M 



Br(/j, ->• ey) V"W r ^ ° 1§ ' ^ 
for both linear and inverse seesaw and Br(r — > fxy) / Br(r — > ey). 

7.2 Supersymmetric high-scale seesaw and flavour symmetry 
In the presence of supersymmetry, the lepton flavour violation observed in neutrino oscillations induces 



decays like /i — > e y, flavour violating tau decays as well as nuclear /i — e conversion (Fig. 14 1 
through the exchange of supersymmetric leptons, as discussed for example, in 1 165 - [167 ] and 1 168| | 169| 



Instead of considering a generic "flavour-blind"supersymmetric high-scale seesaw scheme here we con- 
sider, as an example, the BMV model already introduced in Sec. [3] A detailed numerical analysis of lepton 
flavour violation rates has been performed in Ref. | |170| , using constraints from neutrino oscillation data 
and confronting with lepton flavour violation searches |171| . The allowed parameter space is determined 
by a random search through the multi-dimensional parameter space, keeping all supersymmetric masses 
real and in the range 100 GeV to 1000 GeV. The strongest bounds on lepton flavour violation come from 

£j -> hy. 

The allowed range for the charged slepton parameters is quite restricted. The spectra fall into two 
different groups. The normal hierarchy having two low mass sleptons (~ 150 GeV) and one heavy (above 
~ 500 GeV), and the inverted hierarchy case having two heavy sleptons and one light. In both cases at least 
one slepton mass lies below about 200 GeV, detectable at the LHC. Most points fall into the case of normal 
hierarchy, which often corresponds to a normal hierarchy for the neutrinos as well. The typical case has 
one small and two large mixing angles. Evidently the small mixing angle is needed to suppress the decay 
/i — s- ey. Also the degeneracy of two of the sleptons helps to minimize the LFV As a rule of thumb there 
is at least one pair of sleptons with a mass splitting of less than 40 GeV. 

An important outcome of this study is the prediction for the charged lepton decays £j — > £jj. As seen in 



Fig 16 a lower bound of 10 9 for BR(r — > /17) is found. The is within reach of BaBar and Belle searches. 
Also, BR(/i — > ey) is constrained to be larger than about 10~ 15 and therefore stands good chance of being 



observed in the future 1 156]. The value of tan(/3), also plotted in Fig 16 is constrained to be small. For 



□f r • tr 10 a 



10 r 



10 



tan |1 tan li tan p 

Fig. 16 The predictions for the branching ratios for the processes li — > £jj as a function of tan(/3). 



large tan(/3) the RGE effect destroys the agreement with the solar data. The numerical value of 5 T can not 
be much bigger than the solar mass scale. A rough estimate gives S T ~ 5 x 10~ 4 , corresponding to the 
bound tan(/3) < 10. This agrees with the precise bound found in Fig 16 For small values of tan(/3) the 
threshold corrections dominate and the strongest constraint comes from the bound on BR(fi — > ey). 
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7.3 Low-scale seesaw and lepton flavour violation 

Here we consider the rates for the fi~ —> e~j decay in the framework of the supersymmetric inverse 
seesaw model 1 163 164 1. Fig.[l7]displays the dependence of the branching ratios for fx~ — e~ conversion 
in Ti (left) and \i~ — > e~7 (right) with the small neutrino mixing angle #13, for different values of #12 
(black curve: 9%2 best fit value, blue bands denote 2<r, 3a, 4a confidence intervals for the solar mixing 
angle (9 12 ). The inverse seesaw parameters are given by: M = 1 TeV and fi = 30 eV. The light neutrino 
parameters used are from (2), except for 6*13 which is varied as shown in the plots. The vertical lines 
indicate the sin 2 6*13 values indicated by recent experiments, as well as by the global fit in Ref. fl2|. T hese 
lepton flavour violation rates may be testable in the new generation of upcoming experiments 1 155 156 1. 
For large M the estimates recover those of the standard supersymmetric seesaw. 
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Fig. 17 LFV branching ratios in the supersymmetric inverse seesaw model of neutrino mass (see text). 




7.4 Neutral heavy lepton versus supersymmetric exchange 

The novel feature present in low-scale models and not in the minimal seesaw is the possibility of enhancing 
Br(fi —> e-f) and other tau decays with lepton flavour violation due to neutral heavy lepton (right-handed 
neutrino) versus supersymmetric lepton exchange. In the case where M is low, around TeV or so this 
happens even in the absence of supersymmetry. In this region of parameters the model also gives rise 



to large estimates for the nuclear pT — e~ conversion, depicted in Fig. 14 The latter fall within the 
sensitivity of future experiments | |155) . Note that large lepton flavour violation rates are possible even in 
the massless neutrino limit. The allowed lepton flavour and CP violation rates are, in fact, unsuppressed by 



the smallness of neutrino masses [[152 -154 160 161 1. Finally, for low enough M the corresponding heavy 



leptons could be searched directly at particle accelerators such as LEP already did 1 172 173 1. Prospects of 
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LHC detection are less clear, though they are good in the presence of an extended gauge boson "portal", 
such as right-handed gauge bosons 1 158 159) . 



8 Collider tests: probing neutrino flavour mixing at the LHC 

We now turn to the case of low-scale models of neutrino mass. As an example we consider the case of 
models where supersymmetry is the origin of neutrino mass [43,48], considered in Sec. 2.4 A general 
feature of these models is that the lightest supersymmetric particle (LSP) is unstable, since it is not pro- 
tected by any symmetry. In order to reproduce the masses required by current neutrino oscillation data, the 



LSP is typically expected to decay inside the detector, leaving a displaced vertex 1 174 ■ 177) as seen in the 
left panel in Fig. 19 More strikingly, its decay properties correlate with neutrino mixing angles, as seen 



in the right panel. For example, if the LSP is the lightest neutralino, it is expected to have the same decay 



rate into muons and taus, since the observed atmospheric angle, is relatively close to 7r/4 1 178 - 1 80 1 . This 



Decay length without boost (mm) 





500 1000 1500 2000 El 

m„ (GeV) 

Fig. 19 LSP decays typically produce displaced vertices inside the detector, their semi-leptonic decay branching 
ratios correlate well with the atmospheric mixing angle [176 178], whose low energy determination is illustrated by 
the vertical band |12| . 

opens the tantalizing possibility of testing neutrino mixing at high energy accelerators, like the LHC and 
the "International Linear Collider" (ILC) and constitutes a smoking gun signature of this proposal that for 
sure will be tested. This possibility also illustrates the complementarity of accelerator and non-accelerator 
approaches in elementary particle physics. Before closing this discussion we mention a recent attempt to 
introduce a flavour symmetry to the bilinear R-parity violation scheme [ 181 1. 



9 Implications of a "large" reactor angle 



Recently reactor experiments Double CHOOZ [6|, Daya Bay [7| and RENO (HI have published 
0.092 ± 0.016(stat) ± 0.005(syst) at 5.2ct (DayaBay) 



sin 2 26 



13 



sin 2 26»i 3 = 0.113 ± 0.013(stat) ± 0.019(syst) at4.9cr (RENO) 



(56) 
(57) 



with similar results recently presented at the Neutrino 2012 conference in Kyoto. 

Here we argue that the TBM ansatz can still be taken as a good first order approximation. As we 
discussed in the introduction, in order to have the TBM mixing pattern we need to break separately our 
flavour group (for instance A4) into Z3 in the charged sector and into Z2 in the neutrino sector. Therefore 
the flavour group is completely broken. Since the flavour symmetry leading to the TBM ansatz is in general 
broken we expect deviations from TBM which, in particular, could generate a nonzero reactor angle. 

However many models having TBM at leading order are ruled out because of the recent results which 
indicate that sin #13 ~ Ac where Ac « 0.2. In fact in general we expect that next-to-leading order terms 
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give corrections of the same order Sg to the three angles #13, 6*12 and 6*23. Assuming Sg ~ Ac we have 

sin 2 26»i3 = 0.087 for 5g = 0.15 (58) 
sin 2 26*13 = 0.152 for 5g = 0.20 (59) 

close to the best fits in ( |56] > and p7j ). However the deviations of the solar mixing from its trimaximal values 
sin 2 6*j r 2 BM = 1/3 will be too large if we take Sg = 0.15 ~ Ac, namely 

sin 2 (6*f 2 BM + <5 e ) = 0.48 (0.38 @3<r) (60) 
sin 2 (0f 2 BM -Sg) = 0.20 (0.27 @3<r). (61) 

While this poses no problem for the BMV model, which does not predict the solar angle JT3), it in principle 
rules out most TBM schemes. Indeed, most extensions of TBM models which allow for a large reactor 
angle also predict a deviation of the atmospheric, see for instance |182||183) , and/or the, by now well- 
measured, solar mixing angle from their TBM values. Therefore one of the most relevant theoretical and 
experimental questions is to evaluate the extent to which solar and atmospheric mixing angles deviate from 
their TBM values. 

Still, not all TBM models proposed in the past are excluded, for example in the model of Ref. ]143| , 
based on A4,, large reactor angle 6*13 ~ Ac has been obtained with deviation of 9f 2 BM of order of A c 
in agreement with data. There are other examples in the literature of models where such deviations are 
negligible, despite the relatively large reactor angle value, see for instance | |184f[188) . 

Many alternative ansatze have been considered in order to circumvent this problem. An interesting 
possibility is that the leading order neutrino mass matrix is not diagonalized by the TBM ansatz, but rather 
by the bi-maximal one (where both solar and atmospheric mixing angles are maximal from the start) [ 1 89 1 



or simply bi-large [ 190], or by the golden ratio scheme [ 188 , 191]. Clearly a "large" reactor angle will not 
only act as a "portal" to a new world of CP violation in the lepton sector, but may also shed light into the 
flavour problem, one of the deepest puzzles to our current theories of matter. 
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